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Abstract 



Dark matter candidates arising in models of particle physics incorporating 
weak scale supersymmetry may produce detectable signals through their anni- 
hilation into neutrinos, photons, or positrons. A large number of relevant ex- 
periments are planned or underway. The 'logically possible' parameter space 
is unwieldy. By working in the framework of minimal supergravity, we can 
survey the implications of the experiments for each other, as well as for di- 
rect searches, collider searches, low-energy experiments, and naturalness in 
a transparent fashion. We find that a wide variety of experiments provide 
interesting probes. Particularly promising signals arise in the mixed gaugino- 
Higgsino region. This region is favored by low-energy particle physics con- 
straints and arises naturally from minimal supergravity due to the focus point 
mechanism. Indirect dark matter searches and traditional particle searches 
are highly complementary. In cosmologically preferred models, if there are 
charged superpartners with masses below 250 GeV, then some signature of 
supersymmetry must appear before the LHC begins operation. 
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I. INTRODUCTION 



It is now well established that luminous matter makes up only a small fraction of the 
mass of the observed universe. The evidence for dark matter is both astrophysical and cos- 
mological Such evidence requires only that dark matter is gravitationally interacting. 
However, additional constraints, especially the success of light-element cosmonucleosynthe- 
sis calculations, strongly disfavor the possibility that dark matter is composed solely of 
baryons 0, and so some form of matter foreign to our everyday world is required. The 
dark matter problem is therefore also an important problem for particle physics, as particle 
physics both suggests promising possibilities and imposes stringent constraints. 

Neutralinos are well-motivated candidates to provide much or all of the non-baryonic 
dark matter. An effectively stable particle is a generic component of models with weak- 
scale supersymmetry. This particle is the lightest supersymmetric particle (LSP), and is 
typically the neutralino, a mixture of the superpartners of Higgs and electroweak gauge 
bosons. Particle physics considerations alone require the neutralino to be electrically neutral, 
effectively stable (assuming i?-parity conservation, which is also motivated by the need 
to forbid too-rapid proton decay), and weakly interacting, with mass of order 100 GeV 
(required, as we shall quantify below, if supersymmetry naturally protects the electroweak 
scale from large radiative corrections). Remarkably, these properties are consistent with the 
possibility that the thermal relic density of neutralinos makes up most of the missing mass 
of the universe . 

Unfortunately, these properties also guarantee that neutralinos are practically impossible 
to observe in collider experiments directly. They pass through collider detectors without 
interacting. Existing bounds on neutralinos therefore rely on model-dependent correlations 
between their properties and those of other supersymmetric particles. At present, in minimal 
supergravity, the LEP experiments constrain ^ 40 GeV In the next several years. 
Run II of the Tevatron at Fermilab and, eventually, the Large Hadron Collider at CERN 
will provide more powerful collider probes. 

If neutralinos make up a significant portion of the halo dark matter, many additional av- 
enues for their detection open up. They may deposit energy as they scatter off nuclei in detec- 
tors. We have investigated the prospects for direct detection in a companion article , where 
we emphasized the importance and promise of a mixed gaugino-Higgsino regime, previously 
neglected. Here we will study the possibility of detecting neutralinos indirectly by looking for 
evidence of their annihilation |Q . In the next five years, an astounding array of experiments 
will be sensitive to the various potential neutralino annihilation products. These include 
under-ice and underwater neutrino telescopes (AMANDA, NESTOR, ANTARES), atmo- 
spheric Cherenkov telescopes (STACEE, CELESTE, ARGO-YBJ, MAGIC, HESS, CAN- 
GAROO, VERITAS), space-based 7 ray detectors (AGILE, AMS/7, GLAST), and anti- 
matter/anti-particle experiments (PAMELA, AMS). In many cases, these experiments will 
improve current sensitivities by several orders of magnitude. 

In this paper we evaluate the prospects for neutralino dark matter discovery through 
indirect detection. The neutralino signals depend on many unknown parameters. At the 
same time, an abundance of theoretical and experimental information from particle physics 
can be brought to bear. The implications of traditional particle physics searches for dark 
matter searches, and vice versa, are already significant, and promise to become much stronger 
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over the next few years. One of our main conclusions is that in a class of particle physics 
models favored by current particle physics constraints, astrophysical signals are especially 
enhanced. 

Previous discussions of indirect neutralino detection fall rather sharply into two schools. 
Several previous works are based on specific high energy models , including a number in 
the framework of minimal supergravity [p!0| -[T^ . As we will recall below, minimal supergrav- 
ity incorporates several desirable features, including the radiative breaking of electroweak 



symmetry and the possibility, suggested by gauge coupling unification |T5[, of perturba- 
tive extrapolation to large energy scales. Previous studies in minimal supergravity have 
concluded that only Bino-like dark matter is allowed by particle physics constraints. Such 
dark matter necessarily implies highly suppressed dark matter signals, as we will see. These 
studies, and their somber conclusions, have been criticized as products of overly restrictive 
particle physics assumptions 0]. Recently we have argued more specifically that, even in 
minimal models, these studies failed to examine a very well-motivated regime of parameters, 
and that for this reason their conclusions are overly pessimistic 0. 

At the other extreme, several studies scan over a large set of weak-scale supersymmetry 
parameters and consider values for these parameters as large as 50 TeV. (See, for example, 
Refs. []TB|-^].) These studies, and others, bring a high level of sophistication to the evaluation 
of astrophysical effects on dark matter signals. In this regard, we will have nothing to add, 
but we will incorporate many of the most accurate recent results in our study. 

From a particle physics perspective, this second group of studies is impressively general, 
but this generality is achieved at a cost. For example, extrapolating a given set of weak- 
scale parameters to higher scales, one may encounter such diseases as Landau poles or 
charge- or color-breaking minima. Models of this sort do not arise within a reasonable high 
energy framework. In addition, there is the practical difficulty that the proliferation of 
free parameters implies that results cannot be presented in a systematic, yet transparent, 
fashion. Typically they are displayed as scatter plots after scanning over all parameters. It 
is nearly impossible, from such plots, to determine the dependence of the signal rates on the 
underlying physical parameters. Dedicated correlation plots have been used to highlight a 
few of the relations between dark matter detection experiments, but even the most general 
implications of these experiments for collider searches, electric dipole moments, anomalous 
magnetic moments, proton decay, flavor violation, and other searches for supersymmetry 
are very hard to discern. Finally, it can be expected that supersymmetry parameters, such 
as the /i and gaugino mass parameters, of order 50 TeV will require fine-tuning of the order 
of 1 part in 10^ to produce the observed electroweak scale. (While it is impossible to speak 
of naturalness without first specifying a mechanism of electroweak symmetry breaking, at 
present all concrete models display this feature.) Such large fine-tuning destroys one of 
the main motivations for considering supersymmetric extensions of the standard model in 
the first place. Lacking both a systematic framework and a systematic presentation, it is 
impossible to see how the expectations narrow when some naturalness criterion is imposed. 

It may be impossible to satisfy both schools simultaneously. However, recent results 



suggest an appealing compromise. As has been emphasized in Refs. [pi| -|23[|, a Bino-like 
LSP is not a robust prediction of minimal supergravity. We find that both Bino-like and 
mixed gaugino-Higgsino dark matter is possible. So simply by considering all of minimal 
supergravity parameter space carefully, as we will do here, we remove the most egregious 
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form of model dependence. At the same time, by staying within the confines of minimal 
supergravity we will be able to present results in an organized and comprehensive manner, 
so that correlations with all other supersymmetric signals are easily determined. As the 
experiments discussed below report results, it will be ever more interesting to see what 
models are being excluded or favored. The framework discussed here makes this possible. 

Inclusion of the new gaugino-Higgsino LSP models in minimal supergravity is not just 
a formality. The region with mixed gaugino-Higgsino LSPs is now known to be robust 



and natural, given an objective definition of naturalness [21-23]. It yields cosmologically 



interesting relic densities M, and is even favored by low energy constraints such as proton 



decay and electric dipole moments EM. As we will see, all indirect detection signals are 



enhanced in this region. This lends increased interest to indirect dark matter searches, since 
large — possibly spectacular — rates are predicted within an attractive and simple high 
energy framework. 

In the following section, we review a few essential results concerning neutralino dark 
matter in minimal supergravity with an emphasis on the new gaugino-Higgsino LSP region. 
We then consider each of three promising signals in the following sections: upward-going 
muons from neutrinos in Sec. |ni| , photons in Sec. and positrons in Sec. 0. In Sec. ^ 
we compare these to direct dark matter and traditional particle physics searches, and in 
Sec. |V11| we summarize our results. 



II. NEUTRALINO DARK MATTER AND ITS ANNIHILATION 

The lightest neutralino is the LSP in many supersymmetric models. Assuming i?-parity 
conservation to prevent too-rapid proton decay through dimension-four operators, the LSP 
is effectively stable, and the neutralino is then an excellent candidate for cold dark matter. 

The signals of neutralino dark matter are determined in large part by their composition. 
Neutralinos are mixtures of the superpartners of Higgs and electroweak gauge bosons. After 
electroweak symmetry breaking, these gauge eigenstates mix through the tree-level mass 
matrix 

/ Ml —TTiz cos P sin 9w mz sin [3 sin 6w \ 

M2 mz cos P cos 9w —mz sin (3 cos 9w 

—mz cos P sin 9w f^z cos (3 cos 9w — yU 

\ mz sm. [3 sin 9 w — mz sin P cos 9w — / 

in the basis {—iB, —iW^, H^, Hd)- The weak scale parameters entering this mass matrix are 
the Bino, Wino, and Higgsino mass parameters Mi, M2, and /i, and the ratio of Higgs scalar 
vacuum expectation values tan P = / (H^) ■ The lightest neutralino mass eigenstate is 
then determined by these parameters to be some mixture 

X = xl = a^{-iB) + a^i-iW^) + a-^Hl + 04^,° . (2) 

We define the LSP gaugino fraction to be 

= |aip + |a2p . (3) 
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In the following, we refer to neutralinos with 0.9 < as gaugino-like, 0.1 < < 0.9 as 
mixed gaugino-Higgsino, and R^ < 0.1 as Higgsino-like. 

The preceding discussion is model-independent, assuming only minimal field content. 
However, the minimal supersymmetric standard model is undoubtedly a low-energy effective 
theory of a more fundamental theory defined at some higher scale, such as the grand unified 
theory (GUT) or string scale. A simple realization of this idea is the framework of minimal 
supergravity, which is fully specified by the five parameters (four continuous, one binary) 

mo, Ml/2, v4o, tan/5, sgn(/i) . (4) 

Here, mo, M1/2, and Aq are the universal scalar mass, gaugino mass, and trilinear scalar 
coupling. They are assumed to arise through supersymmetry breaking in a hidden sector at 
the GUT scale Mqut — 2 x 10^^ GeV. The hidden-sector parameters then determine all the 
couplings and masses of the weak scale Lagrangian through renormalization group evolu- 
tion. In particular, electroweak symmetry is broken radiatively by the effects of a large top 
quark Yukawa coupling, and the electroweak scale is determined in terms of supersymmetry 
parameters through 



1 2 _ m%^-ml jan^(3 
2""^" tan^/? 



ml= "\ 2. : -/^ • (5) 



Equation (^ receives corrections beyond tree-level; in our work we include all one-loop effects 
in the Higgs potential We also use two-loop renormalization group equations with 



one-loop threshold corrections p5| , |27| and calculate all superpartner masses to one-loop p5| 



All of the qualitative features to be described below are present, however, also for one- loop 
renormalization group equations. 

Minimal supergravity is, of course, by no means the most general allowed framework. It 
is worth noting, however, that the assumptions most relevant for dark matter, namely, the 
universality of gaugino and scalar masses, are motivated not only by their simplicity, but also 
by concrete experimental facts. The case for gaugino mass unification is especially powerful. 
In the minimal supersymmetric standard model, the three gauge couplings, when evolved 
to high scales, meet with high precision at a point, the GUT scale If the standard 

model is unified in a grand unified gauge theory, one typically obtains also gaugino mass 
unification. The unification of couplings calculation also distinguishes Mqut — 2 x 10^^ GeV 
as the natural scale for a more fundamental framework. 

The unification of scalar masses is motivated by a similar, although more speculative, 
argument. Consider the mass parameter m^^, which, from Eq. (^, plays the critical role 
in determining the weak scale for all moderate and large values of tan/? (tan/3 ^ 5). For 
theories with a universal scalar mass, it is a remarkable fact that the renormalization group 
trajectories of m^^ for various initial conditions mo, when evolved to low scales, meet with 
high precision at a point, the weak scale pT] , p2| . This focusing, which requires that the top 



quark mass be within ~ 5 GeV of its measured value, implies that the electroweak potential 
is highly insensitive to mo. The longstanding problems of supersymmetry with respect to 
CP violation, proton decay, etc. can therefore be ameliorated without fine-tuning, simply 
by assuming large scalar masses. Although the focusing property holds more generally, 
its minimal and most concrete realization is in theories with a universal scalar mass. The 
assumption of a universal (and large) scalar mass is therefore motivated by the fact that it 
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provides a simple and elegant solution to several well-known plienomenological problems of 
weak scale supersymmetry. 

Given these motivations for minimal supergravity, we now consider their implications for 
neutralino dark matter. Gaugino mass unification implies Mi ~ M2/2 ~ 0.4Mi/2- Dark 
matter is therefore never Wino-like,Q and in fact, throughout parameter space, |a2p < 0.07. 
Additional insights follow from re-writing Eq. (^) in terms of GUT scale parameters. For 
Aq = and tan (3 = 10, for example. 



1 

-m 



2 



2 ~ -0.02 ml + 0.7 Ml^ - /i^ , (6) 



^ ~ u.uz, //tg -r u. I ivii/2 



where the numerical coefficients of the first two terms vary fractionally by (9(10%) in the 
(mo, Ml/2 ) plane |]2l| , ^ . The coefficient of ml is highly suppressed |3^j3l|. This is another 



expression of the focusing behavior discussed above, and implies that multi-TeV values of 
mo do not involve significant large fine-tuning. The coefficient is also negative. For fixed 
Ml/2, as mo increases, decreases, and the LSP becomes increasingly Higgsino-like. This 
is important, because even a 10% Higgsino admixture drastically affects the phenomenology. 
In Figs. |l] and ^ we show the LSP mass and gaugino fraction in the (mo, M1/2) plane. For 
large mo ^ 1 TeV, we find a region, previously ignored, where the LSP has a significant 
Higgsino component. The green shaded regions are excluded by the requirement that the 
LSP be neutral (top left) and by the chargino mass limit of 95 GeV (bottom and right). 

In Figs, m and ^ we also indicate the regions of cosmologically interesting relic densities; 
see Ref. for details. A generous range is 0.025 < Q^h"^ < 1, where the lower bound is 
the requirement that neutralino dark matter explain galactic rotation curves and the upper 
bound follows from the lifetime of the universe. Above this shaded region, Qh"^ > 1; below, 
Qh'^ < 0.025. The range 0.1 < f^x^^ — ^-^ is most preferred by current limits. Our relic 
density calculation is not trustworthy in the black region, where there is a Higgs scalar 
resonance, and very near to the left and right borders of the excluded region, where co- 



annihilation is important ||3^-p5| . In the bulk of parameter space, however, these effects are 
negligible. For all tan/5, cosmologically interesting densities are possible for mo ^ 1 TeV. 
For tan/3 = 10, the cosmologically preferred region contains gaugino-Higgsino dark matter. 

In contrast to mo, the parameters M1/2 and /i enter Eq. (|^) with 0{1) coefficients. Natu- 
ralness therefore requires that the LSP mass (and, in fact, the masses of all four neutralinos 
and both charginos) should not be too far above the electroweak scale. While in principle 
it is possible that in some fundamental framework M1/2 and /i are correlated precisely in a 
way that allows both parameters to be large without fine-tuning (a possibility considered in 
Ref. p6[), no such framework has been found to date. Barring such a possibility, extreme 
values such as Mi/2,/i ~ 50 TeV require a fine-tuning of 1 part in 10^ and destroy one of 
the prime motivations for weak-scale supersymmetry. We therefore regard such large values 
as highly disfavored, and we will focus on neutralino masses of order 100 GeV. 



Wino-like LSPs exist in other frameworks, but typically they annihilate far too quickly to be 
cosmologically relevant Interesting relic densities are possible, however, if there is some 

mechanism of late production ||2^ . 



6 




FIG. 1. Contours of constant LSP mass in GeV in the (mo,M^/2) plane for Aq = 0, fj, > 0, 
rrit = 174 GeV, and two representative values of tan /3. The green shaded regions are excluded by 
the requirement that the LSP be neutral (left) and by the chargino mass hmit of 95 GeV (bottom 
and right). We have also delineated the regions with potentially interesting values of the LSP relic 
abundance: 0.025 < ^x^'^ — ^ (yellow) and 0.1 < i^^/i^ < 0.3 (light blue). In the black region, 
|2m^ — nifil < 5 GeV, and neutrahno annihilation is enhanced by a Higgs resonance. 




mo (GeV) mo (GeV) 

FIG. 2. Contours of constant gaugino fraction in percent, for the same values of the pa- 
rameters as in Fig. |l[ 

Neutralinos annihilate through a variety of channels. The three leading processes are 
shown in Fig. ^. (Note that co-annihilation, while potentially important in determining 
relic densities in the early universe, is negligible now.) Annihilation into gauge bosons is of 
particular importance, as these processes lead to more energetic and striking signals. The 
WW cross section relies on Wxxf interactions. The only such couplings allowed by gauge 
invariance are WH^H"^ and WW^W^. However, as noted above, gaugino mass unification 
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FIG. 3. Three leading neutralino annihilation channels. 
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FIG. 4. Contours of constant av in pb for (a) XX ~^ WW and (b) XX ~^ We fix Aq = 0, 
0^ mt = 174 GeV, and tan/3 = 10. 



implies that the Wino content of the LSP is always negligible. A large WW cross section 
is therefore possible only when the LSP has a significant Higgsino component. The same 
conclusion holds for the ZZ process, where the Zxx^ interaction is possible only through 
ZH^H^ couplings. In Fig. ^, we see that the annihilation cross sections for xx WW and 
XX ~^ 8^re indeed highly suppressed in regions with Bino-like LSPs, but are enhanced 
by three to four orders of magnitude in regions with mixed gaugino-Higgsino dark matter. 
As we will see, this region, favored by low energy constraints, will be the most promising 
for all indirect signals. 

Before closing this section, we note several features of Figs. || and |^ that will also apply to 
many of the following figures. Unless otherwise noted, we present results for Aq = 0, ft > 0, 
rrit = 174 GeV, and representative values of tan (3 as indicated. Aq governs the left-right 
mixing of scalars, and does not enter the neutralino sector. It is therefore largely irrelevant, 
especially in the regions of parameter space with observable signals — where, as we will see, 
the scalars are heavy and decoupled.0 (Besides, the most important trilinear coupling. At, 



^For ^0 = 0, the 1999 LEP bound on the mass of the lightest CP-even Higgs boson rrih > 
107.7 GeV barely constrains the parameter space shown in our figures, while the 2000 Hmit ruh > 
113.3 GeV excludes the region with niQ ^ 1 TeV and M1/2 ~ 300 GeV. However, the Higgs search is 
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has a weak-scale fixed point, and so is only weakly sensitive to Aq.) Our dark matter results 
are rather insensitive to the sign of fi, but the choice /i > is motivated by the constraint 
from B — >• (see, e.g. [0). Finally, perturbativity of Yukawa couplings limits tan/? to 
the range 1 < tan/5 ^ 60. Low values of tan/? < 3 are now being excluded by the LEP 
Higgs search. In the remaining interval, models with moderate and high values may have 
qualitatively different behavior, as processes proportional to down-type Yukawa couplings 
are enhanced by tan /?. We therefore typically present results for two representative values, 
one in each range. Plots of many other quantities, including all physical slepton, squark, 
and Higgs masses in the (mo,Mi/2) plane, including the high mo region, may be found in 
Refs. |21^. 

Finally, in the following three sections, we will not include the effects of variations in 
thermal relic density in our signal rates, but rather assume, for concreteness, a fixed local 
neutralino density. The results are then more transparent, and are applicable to general 
scenarios, such as those in which a late source of neutralino production is present. Of course, 
in the simplest scenario, models with f^^^^ ^ -'- excluded, and those with an under- 
abundance of neutralinos are disfavored and imply suppressed or negligible dark matter 



signals. This should be kept in mind in the following sections. In Sees. and [VI I| , we will 
combine all these considerations, and focus on the most preferred regions. 



III. NEUTRINOS 

When neutralinos pass through astrophysical objects, they may be slowed below escape 
velocity by elastic scattering. Once captured, they then settle to the center, where their 
densities and annihilation rates are greatly enhanced. While most of their annihilation 
products are immediately absorbed, neutrinos are not. High energy neutrinos from the 



cores of the Earth p8|- ^] and Sun p0| , ^2| -|47|| are therefore promising signals for indirect 
dark matter detection. 

The formalism for calculating neutrino fluxes from dark matter annihilation is well de- 
veloped. (See Ref. |^ for a review.) The neutrino flux depends first and foremost on 



the neutralino density, which is governed by the competing processes of gravitational cap- 
ture and neutralino annihilation. If N is the number of neutralinos in the Earth or Sun, 
N = C — AN"^, where C is the capture rate and A is the total annihilation cross section 
times relative velocity per volume. The present neutralino annihilation rate is then 

= ^AN' = ^Ctfuih^VCAt^) , (7) 

where ~ 4.5 Gyr is the age of the solar system. 

Captured neutralinos then annihilate through the processes of Fig. (H). As XX ~^ ff 
is helicity-suppressed, neutrinos are produced only in the decays of primary annihilation 



highly sensitive to the Aq parameter. For other values of the current constraints may be evaded 
with negligible impact on our predictions for indirect dark matter detection. For this reason, we 
do not include the Higgs mass constraints in the following analysis. 
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products. Typical neutrino energies are then Ei, ~ -^m^ to with the most energetic 

spectra resulting from WW, ZZ, and, to a lesser extent, rf . After propagating to the 
Earth's surface, neutrinos are detected through their charged-current interactions. The 
most promising signal is from upward-going muon neutrinos that convert to muons in the 
surrounding rock, water, or ice, producing through-going muons in detectors. The detection 
rate for such neutrinos is greatly enhanced for high energy neutrinos, as both the charged- 
current cross section and the muon range are proportional to E,^. 

The calculation of muon fluxes from neutralino annihilation in the Earth and Sun is on 
reasonably firm footing, as it depends only on the local dark matter density and is insensitive 
to details of halo modeling.^ Nevertheless, the calculation is involved, primarily as a result of 
complications in evaluating capture rates l^l]-^ and, in the case of the Sun, in propagating 
the neutrinos from the core to the surface p0|,|54|,p5| . Here we make use of the procedure 



of Refs. Il48|j56[| . For other analyses, see Refs. |l57|- ^ , |l^ , as well as those motivated by the 



Tevatron e^e~77 event []60 , |6T| an d by the DAMA annual modulation signal |62,K 

The 'filling parameters' VCAtQ for the Earth and Sun are given in Figs. || and ^ For the 
Sun, y/CAtQ ^ 1 for all supersymmetry parameters. The neutralino density has therefore 
reached equilibrium, and the annihilation rate is at full strength, with ^ C/2. For the 
Earth, however, typically VCAIq -C 1, and the annihilation rate is F^ ^ ^C'^At^ and far 
from maximal. As we will see, this plays an important part in reducing the Earth's signal 
below the Sun's. 

The other major ingredient in the muon flux computation is the estimate of the neutralino 
capture rate C, which is shown in Figs. ^ and |^ for the Earth and Sun, respectively. The 
elemental compositions of the Earth and Sun are given in Ref. A quick comparison of 
Figs. 1^ and || reveals that, not surprisingly, a large astrophysical body like the Sun is much 
more efficient in trapping neutralinos. The Earth's capture rate is, however, enhanced by the 
iron resonance for very light neutralino masses ~ 50 GeV. The tan j3 dependence is also 
noteworthy. The capture rate in the Earth is determined primarily by the spin-independent 
elastic scattering cross section for XQ ~^ XQ through s-channel squarks and t-channel Higgs 
boson exchange. All amplitudes require chirality flips, either through Higgs interactions, 
squark mass insertions, or quark mass insertions. The first two are proportional to tan/3 
and therefore dominate for moderate and large tan (3, leading to C ~ tan^ (3. In contrast, for 
the Sun, the tan/3 dependence is minimal; the dominant contribution is from axial- vector 
scattering off Hydrogen, which is largely independent of tan f3. 

Muon flux rates from the Earth and Sun are presented in Figs. |^ and Consistent 
with previous studies (see, e.g. [|14|)) we flnd that the flux rate is indeed small in regions of 



parameter space with mg < 1 TeV and Bino-like LSPs. However, for uiq > 1 TeV, in the 
region where > mw and the dark matter is a gaugino-Higgsino mixture, the fluxes are 
greatly enhanced. Here, annihilation to gauge bosons is unsuppressed, resulting in a hard 
neutrino spectrum and large muon fluxes. In this region, the rates from the Sun are large 



It has been suggested that muon fluxes may be enhanced, by up to two orders of magnitude, 
due to capture of neutralinos in highly eccentric solar system orbits The magnitude of the 

enhancement depends on details of the neutralino parameters and involves astrophysical issues still 
under debate We have not included it here. 



10 




for all values of tan/3. For the Earth, we see that, despite the close proximity of the Earth's 
center, the muon fluxes are typically suppressed by several orders of magnitude relative to 
those from the Sun. However, reasonably large rates are possible even for the Earth for large 
tan /3, and particularly for very light neutralinos, where the capture rate is enhanced by the 
iron resonance, as discussed above. 



600 
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FIG. 6. As in Fig. |5|, but for the Sun. 



The theoretical predictions of Figs. |^ and |T0| can be compared with the experimental 
sensitivities of ongoing and near future detectors [Q. These experiments, along with their 
more salient characteristics and flux limits (where available), are listed in Table |. The 
flux limits depend on the expected angular dispersion in the signal. This dispersion has 
two possible origins. One is the source: although neutralinos from the Sun are essentially 
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FIG. 7. The capture rate of neutralinos C in in the center of the Earth (for v 
and pQ = 0.3 GeV/cm^). 



270 km/s 



a point source, in the Earth, 98% of neutrahno annihilations occur within a cone of half- 
angle 8.6°y^50 GeV/m^ 1111. The second is the angle ^^ms ~ 13° -^25 GeV/E^ between 
the neutrino and its daughter muon. As < 6'rms is typically the dominant effect. 

The flux limits listed are for half-cone sizes of 15°, corresponding roughly to m-^ ~ 50 GeV. 
Heavier neutralinos will produce more coUimated muons, allowing smaller cone sizes with 
reduced backgrounds. The improved limits for smaller cone sizes may be found in the 
references. 
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FIG. 8. Same as in Fig. but for the Sun. 



Comparing Figs. ^ and |TD| with Table |, we find that present limits do not significantly 
constrain the minimal supergravity parameter space. However, given that the effective area 
of neutrino telescope experiments is expected to increase by 10 to 100 in the next few years. 
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muon fluxes of order 10-100 km^^ yr^^ may be witliin reacli. Sucli sensitivities are typically 
not sufficient to discover Bino-like LSPs, unless they are light and tan (3 is large. But they 
have an excellent opportunity to detect dark matter in the mixed gaugino-Higgsino dark 
matter scenarios, which, as we have emphasized above, are preferred by low energy particle 
physics constraints. 
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FIG. 10. As in Fig. 
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Muon energy thresholds, listed in Table |, are not included in Figs. ^ and |T0[ Since 
the muon detection rate is dominated by high energy muons as noted above, the threshold 
energy is typically not important, especially in the regions where a detectable signal is 
expected. This is not the case for all detectors, however. For example, since muons lose 
0.26 GeV per meter in water and ice, neutrino telescopes requiring track lengths of ~ 100 
m will have thresholds of order ~ 25 GeV. The dependence on threshold energy has been 
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TABLE I. Current and planned neutrino experiments. We list also each experiment's (ex- 
pected) start date, physical dimensions (or approximate effective area), muon threshold energy 



E^}"' in GeV, and 90% CL flux limits for the Earth 
angle ~ 15° when available. 



and Sun $® in km ^ yr 



for half-cone 



Experiment 


Type 


Date 


Dimensions 


J?thr 


$® 

A* 


$0 



Baksan |65] 



Kamiokande [36 
MACRO 11 
Super-Kamiokande |68| 
Baikal NT-96 [16911 



AMANDA B-10 |70[ 
Baikal NT-200 
AMANDA II |7l| 
NESTOR§ @ 
ANTARES [ff 



IceCube |71 



Underground 

Underground 

Underground 

Underground 

Underwater 

Under-ice 

Underwater 

Under-ice 

Underwater 

Underwater 

Under-ice 



1978 
1983 
1989 
1996 
1996 
1997 
1998 
2000 
2000 
2003 
2003-8 



17 X 17 X 11 m^ 
~ 150 m2 
12 X 77 X 9 m^ 
~ 1200 m^ 
~ 1000 m^ 
~ 1000 m^t 
~ 2000 m2 
~ 3 X 10^ m2 
~ 10^ m2t 
~ 2 X 10^ m^t 
~ 10^ m^ 



1 

3 

1.5^ 
1.6 
10 

~ 25 
~ 10 
~ 50 
few 
5-10 



6.6 X 10^ 
10 X 10^ 
3.2 X 10^ 
1.9 X 10^ 

15 X 10^ 
44 X lO^I" 



7.6 X 10^ 
17 X 10^ 
6.5 X 10^ 
5.0 X 10^ 



2 GeV for Sun. "I" Hard spectrum, m^^, = 100 GeV. ^ One tower. 



100 GeV. 



studied in Refs. [0,|l^], where it was found that for threshold energies of E]^ 



~ to 



-m 



the loss of signal is substantial. Low threshold energies in neutrino telescopes are 
clearly very important for dark matter detection. This conclusion is further strengthened 
by considerations of naturalness, which favor low neutralino masses. 



IV. PHOTONS 

High-energy photons provide a unique signal of dark matter annihilation. They point 
back to their source, and their energy distribution is directly measurable, at least in principle. 
For these reasons, given sufficient angular and energy resolution in gamma ray detectors, a 
variety of signals may be considered. 



The photon signal may arise from the galactic center [ffSHfTl, the galactic halo [|78| , |79 [ 



or even from extra-galactic sources [^]. We will consider the galactic center, where large 



enhancements in dark matter density are possible [[19| , p0| . In contrast to the neutrino sig- 
nal considered in Sec. the photon flux is highly sensitive to halo model parameters. 
Fortunately, the problem may be separated into two parts: one containing all halo model 
dependence, and the other all particle physics uncertainties. Given the predicted photon 
fluxes for a reference halo model, the predictions for all other halo models are then easily 
determined. 

The photon energy distribution receives two types of contributions: line and continuum. 



The former results from the loop-mediated processes XX ~^ 77 [eU'EH XX ~^ 1^ ^3 
Because dark matter in the halo is extremely non-relativistic, photons from these processes 
have an energy width of only AE^/E^ ~ 10~^ and are effectively mono-energetic. While 
this signal would be the most spectacular of all possible indirect signals, its rates are, of 
course, suppressed [^. In a model- independent survey, Bergstrom, Ullio, and Buckley [|19 
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have found that the photon hne may be observable for neutrahnos with a large Higgsino 



component, assuming a cuspy halo profile, such as that of Navarro, Frenk, and White [^5 
and telescopes with small angular acceptances ~ 10~^ sr. 

On the other hand, photons may also be produced in the cascade decays of other primary 
annihilation products. In contrast to the line signal, cascade decays produce a large flux 
of photons with a continuum of energies. This signal is far less distinctive and will almost 
certainly require additional confirmation to unambiguously distinguish it from background 
or other exotic sources. Nevertheless, we will focus here on the continuum signal, as it will 
provide the first hint of dark matter from gamma ray astronomy. 

The differential photon flux along a direction that forms an angle ■?/' with respect to the 
direction of the galactic center is 

^ dK 



t A. 

where the sum is over all annihilation channels and p is the neutralino mass density. All of 
the halo model dependence is isolated in the integral, which, following Ref. |]19[, we write in 
the dimensionless form 



"'^''-aJk;^i o.3GeV/cm» j l/^' ■ '^'^ 

The integral is along the line of sight. Assuming a spherical halo, the mass density is given 
by p = p{r), where = l"^ + — 21Rq cosip, and Rq ~ 8.5 kpc is the solar distance to the 
galactic center. 

The photon flux is, of course, maximized for ip = 0, but it must be averaged over the 
field of view. The result is 

$^(Ethr) = 5.6 X 10-^° cm-2 X V HdE^ ( ^] ( ^^^-^]' J{An) Afi , (10) 

■ J Bthr dE \ph J \ rriy^ ) 

where 

J(A^])^-i-/ (11) 

and AQ is the solid angle of the field of view centered on = 0. Eth_r is the lower threshold 
energy; detectors also have upper cutoffs, but these are typically irrelevant, as the energy 
distribution falls steeply with energy. J has been studied for a variety of halo models in 
Ref. |]l9l. For a typical atmospheric Cherenkov telescope (ACT) acceptance of AQ = 10~^ sr, 
the modified isothermal profile described by p(r) oc [1 + (r/a)^] yields 3 ^ J ^ 10^. On 
the other hand, cuspy halos lead to values of J as large as 10^. (Such singular profiles have 
recently been argued to be incompatible with neutralino dark matter, however, based on 
radio emission from neutralino annihilation near the black hole at the galactic center [p6[].) 
We will choose a moderate reference value J(10~^) = 500, which is within the allowed 
ranges of both the modified isothermal and cuspy halos. The factorizability of the photon 
flux implies that our results can be scaled to all other halo models easily. 
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500 1000 1500 2000 500 1000 1500 2000 2500 
mo (GeV) mo (GeV) 

FIG. 11. Photon flux <I>^(i?thr) in cm^'^ from a 10^^ sr cone centered on the galactic center 
for a threshold energy of -Ethr = 1 GeV. We assume halo model parameter J = 500. Results for 
other halo models may be obtained by scaling to the appropriate J (see text). 
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TABLE II. Some of the current and planned 7 ray detector experiments with sensitivity to 
photon energies 10 GeV < i?^ < 300 GeV. We list each experiment's (proposed) start date and 
expected coverage in GeV. The energy ranges are approximate. For experiments constructed 
in stages, the listed threshold energies will not be realized initially. See the references for details. 



Experiment 



Type 



Satellite 
ACT array 
ACT array 
Air shower 
ACT 
Satellite 
ACT array 
Space station 
ACT array 
ACT array 
Satellite 



Date 



Range 



EGRET 
STACEE |8| 
CELESTE @ 
ARGO-YBJ [|l| 
MAGIC i2i 



AGILE |93| 
HESS in 
AMS/7 [H 
CANGAROO III 



VERITAS 



GLAST |98| 



1991-2000 

1998 

1998 

2001 

2001 

2002 

2002 

2003 

2004 

2005 

2005 



0.02-30 
20-300 
20-300 
100-2,000 
10-1000 
0.03-50 
40-5000 
0.3-100 
30-50,000 
50-50,000 
0.1-300 



The particle physics model dependence enters through all the other factors of Eq. (^. 
The energy integral is roughly / dE dN^/dE ~ 0.5 for all i, but the energy distribution de- 
pends significantly on the annihilation channel. The differential gamma ray multiplicity has 
been simulated for light and heavy neutralinos in Refs. and ||19| , respectively. The spec- 
trum for the most important annihilation channels is described well by dN^/ dx = a e~^^ / x^'^, 
where x = E^/m^ and (a, b) = (0.73, 7.76) for WW and ZZ [|T9[, (1.0, 10.7) for bb, (1.1, 15.1) 
for tt, and (0.95, 6.5) for uu. We neglect Higgs boson final states, as they never have branch- 
ing fraction greater than 7%. For the gg final state, we use the light quark distribution. 
Our crude approximation for gluons is relevant only in isolated regions with Bino-like LSPs 
where, as we will see, the signal is unobservable. With the exception of the irrelevant light 
quark distribution, the WW and ZZ distributions produce the most energetic photons. 

The photon flux $^(ii^thr) is given in Figs. 0-|TB| for threshold energies of 1, 10 and 50 
GeV. The maximal rates are found in the region of parameter space with mixed gaugino- 
Higgsino dark matter, and are insensitive to tan/5. Here branching ratios to gauge bosons 
are large, and the photon spectrum hard. In the rest of parameter space, bb is an important 
final state, and $^(i?thr) is enhanced by tan/?. 

In the past, the high energy photon spectrum with 10 GeV ^ E^ ^ 300 GeV has been 
largely unexplored. Ground-based detectors, such as the Whipple 10m telescope, have large 
effective areas, but have traditionally been limited to energies above ~ 300 GeV. Space- 
based detectors, such as EGRET, have been sensitive to photon energies up to ~ 20 GeV, 
but are limited above this energy by their small effective area. There has therefore been an 
unexplored gap at intermediate energies, which happens to overlap substantially with the 
range of energies most favored by supersymmetric dark matter. 

The experimental situation is changing rapidly, however. Currently, two heliostat arrays, 
STACEE and CELESTE, are running with sensitivity in the range 20 GeV ^ E^ ^ 300 GeV, 
and many more experiments with greatly improved sensitivity are expected in the next few 
years. Upcoming experiments with sensitivity to 7 rays with 10 GeV ^ E^ ^ 300 GeV are 
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Ethr (GeV) 

FIG. 14. Integral photon 

fluxes ^^{Ethj.) as a function of the threshold energy E'thr for Aq = 0, fj. > 0, mt = 174 GeV, 
and halo parameter J = 500. The four models have relic density ^x^^ ~ 0.15, and are specified 
by {tanp,mo,Mi/2,mx,Rx) = (10,100,170,61,0.93) (dotted), (10,1600,270,97,0.77) (dashed), 
(10, 2100, 500, 202, 0.88) (dot-dashed), and (50, 1000, 300, 120, 0.96) (solid), where ah masses are in 
GeV. Point source flux sensitivity estimates for several gamma ray detectors are also shown. (Care 
should be taken in comparing these sensitivities to the predicted fluxes — see text.) 



listed in Table ||. 

An important figure of merit for the detection of 7 rays from the galactic center is the 
point source flux sensitivity. A compilation of previous estimates of flux sensitivities is given 
in Fig. m for EGRET, STACEE, CELESTE, ARGO-YBJ, MAGIC, AGILE, HESS H, 
AMS/7 [0, VERITAS [0, and GLAST |8[. The flux sensitivities for the first six exper- 



iments are from Ref. ||93|, and those for the remaining experiments are from the references 
listed. The sensitivity of MAGIC assumes the availability of high quantum efficiency photo- 
sensors. The sensitivity for CANGAROO III is currently being re-evaluated ||100|| . We have 
included it in accord with expectations that it will be comparable to that of HESS. The point 
flux sensitivities are, of course, dependent on the source's location and energy spectrum. 
They are also subject to a variety of other experimental uncertainties and assumptions; see 
the references for details. A typical estimate ||9^ assumes background extrapolated from 
EGRET data and a signal distribution dN^/dE oc E'"^'^. Detector efficiencies and cuts 
are included, and a 5a signal with at least 10 photons is required. 50 hours of observation is 
assumed for telescopes, and one year of an all sky survey for the space-based detectors. The 
arrow for AMS/7 Fig- indicates that a published estimate exists only for E'thr = 1 GeV, 
but flux sensitivity at some level can be expected out to the detector limit of 100 GeV. 

The experimental sensitivities presented clearly cannot be interpreted as future dark 
matter discovery contours. The neutralino signal has a different energy spectrum than 
assumed, and the background in the direction of the galactic center is larger, due to the 
diffuse 7 ray emission from the galactic disk ||101||, which enhances ^/B by a factor of ~ 5 
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102|| . (This last fact implies that for some halo profiles, it may be advantageous to center 
the field of view away from the galactic center. This optimization may significantly reduce 
potential losses in signal significance.) In addition, there are many ambiguities in background 
calibration, and, as noted above, the continuum signal is not sufficiently distinct for a simple 
excess to identify the source as dark matter annihilation. However, the flux sensitivities of 
Fig. |l^ do clearly portend substantial progress in the next few years, and can serve as rough 
indications of what signal levels will be detectable. 

The expected fluxes for four supersymmetry models and J = 500 are also shown in 
Fig. I^. Although there is a large uncertainty from halo model dependence, it is clear 
that detectable signals are possible. At present, EGRET data is not overly constraining, 
although halo models with large J ~ 5000 are within EGRET sensitivity and may even 
explain a flattening of the spectrum. In the future, AMS/7 and AGILE will improve this 
sensitivity, and MAGIC may see excesses for J ~ 500. Finally, GLAST will provide the 
greatest sensitivity of all, probing halo models with J as low as J ~ 50. 

If a significant excess is found in future experiments, its dark matter origin can be tested 
in a variety of ways. Confirmation from other searches for dark matter or supersymmetry 
would be the most satisfying possibility. As we will see in Sec. the neutrino and positron 
signals probe similar models, so a coincidence of various signals is a distinct possibility. 
However, it has also been suggested that an angular distribution of photons that does not 
follow the galactic disk and bulge may be a powerful diagnostic ||19|. Also, as the dark 
matter signal has a shape differing from the background, detailed likelihood fits to the 
photon energy distribution may also be a useful tool, although far beyond the scope of 
this work. It seems clear, in any case, that for reasonable halo models and supersymmetry 
parameters, meaningful 7 ray signals in the next few years are possible, particularly with 
gaugino-Higgsino dark matter. 



V. POSITRONS 



An excess of cosmic anti-particles and anti-matter from dark matter annihilation may be 
detected in space-based or balloon-borne experiments. The positron signal is perhaps the 
most promising ||103| - |107[] . In the past, soft anti-protons with energies ~ 100 MeV have also 
|108| , [L09|] . However, recent work finds larger background than previously 

Anti-deuterium 



been considered 



expected, complicating the identification of an anti-proton signal [110, 111 
has also been suggested clS db possibility [ |112| ]. 

The positron background is most likely to be composed of secondaries produced in the 
interactions of cosmic ray nuclei with interstellar gas, and is expected to fall as ~ 
At energies below 10 GeV, however, this background is subject to large uncertainties from 
the effects of the solar wind [ 106 ,107]. The soft positron spectrum also varies depending 
on the orbit path of the experiment. At high energies, these effects are negligible. In 
addition, positrons lose energy through a variety of processes, and so hard ones must typically 
be produced within a few kpc [|106| , |107|[ . For this reason, the hard spectrum is relatively 
insensitive to variations in the halo profile near the galactic center. The dark matter signal 
is therefore most promising at high energies, where the background is relatively small and 
well understood. 

The differential positron flux is |107[ 
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^ = 4E-.>'a:./<'^^o/,(^;o)G(£„,E), (12) 

where po is the local neutralino mass density, the sum is over all annihilation channels, and 
BI+ is the branching fraction to positrons in channel i. The source function f{Eo) gives the 
initial positron energy distribution from neutralino annihilation. G{Eq, E) is the Green's 
function describing positron propagation in the galaxy ||113|| and contains all the halo model 
dependence. 

For the reasons mentioned above, processes yielding hard positrons are by far the most 
important for dark matter discovery. The 'positron line' signal from e+e~ is helicity- 
suppressed. It may be enhanced, for example, in the case of Bino-like LSPs if selectrons are 
much lighter than all other scalars, but this possibility is highly unmotivated. To an excel- 
lent approximation, then, hard positrons arise from xx ~^ WW, ZZ, followed by the direct 
decay of gauge bosons to positrons. Assuming unpolarized gauge bosons, / is the familiar 
flat distribution with endpoints determined by the gauge boson and neutralino masses. The 
Green's function G has been modeled by Moskalenko and Strong in Ref. |1U7|| in a framework 



that consistently reproduces a wide range of observational data from anti-protons, nuclei, 
electrons, positrons, and photons. 

Combining all of these results, the differential positron flux may be written as 

-4- = 0.027 cm-2 s"^ sr"^ GeV 



X 



po \VlOOGeVV^ aiV 



O.SGeV/cm^J \ J ph ■ (3,, 

where 

1/2 



T.^BUj^"dzg{z,E/m^), (13) 



(3ww,zz = {l -ml^^z/^l) ' (14) 

4 = (l±A)/2, (15) 

B^^ = B{W+ ^ e+z/) = 0.11 , (16) 

5f/ = 2B{Z e+e") = 0.067 , (17) 

and the reduced Green's function is 
g{z,E/m^) = io^ios%E+biog,,E+c _ E/m^) + iQ^^ogl, e +x log,, E+y e{E/m^ - z) , (18) 

where E is in GeV and the (2;-dependent) coefficients a, 6, c and w, x, y are tabulated in 
Ref. ||107|| for different halo profiles. As mentioned above, at high energies, these coefficients 



are fairly independent of the halo model, as high energy positrons originate in our solar 
neighborhood, where all profiles give similar densities. We adopt coefficients corresponding 
to the modified isothermal distribution with halo size 4 kpc. For large m^, the integral 
of Eq. (p!3| ) is insensitive to m^, and so the differential positron flux scales as ~ Xjm^. 
Neutralinos with mass not far above niw are therefore most easily detected. 

In Fig. |T^, we show three sample spectra for supersymmetry models yielding relic abun- 
dances of n^/i^ ^ 0.15. Two background spectra from Ref. ||107|| are also shown. The signal 



rates are signiflcantly suppressed relative to those of Refs. ||105| , |107|| , where the dark matter 
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FIG. 15. The differential positron flux for three of the four sample points in Fig. |lj. The curves 
labeled C and HEMN are background models from Ref. |107]; the dotted curve is our fit to C. 



was assumed to be Higgsino-like. Higgsino-like dark matter is highly disfavored, however, 
as, unless it is unnaturally heavy, it annihilates too strongly to leave interesting relic abun- 
dances. As is evident from Fig. ^, in the allowed minimal supergravity parameter space the 
LSP is far from pure Higgsino-like, particularly in the region with preferred relic density. 
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FIG. 16. The optimal positron energy E'opt in GeV at which the signal to background ratio 
SjB is maximized. 

As positrons result from two-body decay, we expect the signal, and the signal to back- 
ground ratio SjB^ to be peaked near vn^jl. This is evident in the three examples give in 
Fig. 15. In Fig. we plot the optimal energy i?opt at which the signal to background ratio 
is maximized. Our fit to background C is E'^ f.+ / dTldE = 1.16 x lO^^E'^ "^^, where E is 
in GeV. Comparing with Fig. we see that Eopt is indeed approximately itl^/2 through- 
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FIG. 17. The positron signal to background ratio S/B at E'opt- 



TABLE III. Recent and planned detector experiments. We list each experiment's (expected) 
start date, duration, geometrical acceptance in cm^ sr, maximal Ef,+ sensitivity in GeV, and 
(expected) total number of e"*" detected per GeV at = 50 and 100 GeV. 



Experiment 




Type 


Date 


Duration 


Acceptance 




^(50) 


^(100) 


HEAT94/95 


m 


Balloon 


1994/95 


29/26 hr 


495 


50 






CAPRICE94/98 [|TT| 


Balloon 


1994/98 


18/21 hr 


163 


10/30 






PAMELA |T§ 


Satellite 


2002-5 


3 yr 


20 


200 


7 


0.7 


AMS-02 |117 




Space station 


2003-6 


3 yr 


6500 


1000 


2300 


250 



out parameter space. In Fig. |T7|, we plot S/B at i^opt- S/B is substantial only in the 
gaugino-Higgsino region with > mw- 

Figs. |TB| and [T^ imply that the best experimental hope for indirect detection of dark 
matter through positrons is in experiments sensitive to positron energies above ~ 50 GeV. 
In the next two to three years, both PAMELA, a satellite detector, and AMS-02, an exper- 
iment to be placed on the International Space Station, will satisfy this requirement. These 
experiments and other recently completed experiments are listed in Table p| . 

The expected number of positrons per GeV are listed in Table |I| at positron energies 
of 50 and 100 GeV. After integrating over some appropriate energy bin size, we see that 
the expected statistical errors are roughly ~ 10% for PAMELA, and ~ 1% for AMS-02. Of 
course, the signal will also be degraded by systematic errors, particularly in the background 
calculation. It seems likely, however, that the characteristic peaking of the dark matter 
signal over an interval of (9(10 GeV) will be distinctive. In addition, some systematic errors 
may be eliminated by considering the ratio e~^/{e~ + e"*"). While the positron signal is 
typically too small for Bino-like LSPs, an excess of ~ 1% is possible for gaugino-Higgsino 
dark matter. The region of detectable positron signals may be extended, however, if, for 
example, the halo is clumpy, or if the local density is larger than our reference value of 
0.3 GeV/cm^ M] . 
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VI. COMPARISON WITH OTHER SEARCHES 



In the previous three sections, we have examined several indirect signals of neutralino 
dark matter. As emphasized in Sec. |I|, however, supersymmetric dark matter cannot exist 
in isolation, and there are many other avenues for probing supersymmetric models. We 
now discuss several other promising probes and their projected sensitivities, and we then 
compare their reaches. 

Most closely linked to indirect searches are searches for dark matter scattering off nuclei 
in low-background detectors. The DAMA collaboration has reported evidence for an annual 



modulation signal | ]ll(j|| , and the activity in this field will intensify tremendously in the next 
few years. (For a recent review, see, e.g., Ref. ||119|| .) Here we will use our previous results 



to estimate the sensitivities of the direct searches. We choose CDMS (Soudan) | ]120| | and 
CRESST ||121|| as examples of near-term future experiments. Their projected sensitivities 
in neutralino-proton cross section are both of order cXproton ~ 10~^ pb for 50 GeV < < 
500 GeV. More precisely, we parameterize their sensitivities as 



Cproto 



exp <a + b 



100 GeV 



+ c 



100 GeV 



pb 



(19) 



with (a, b, c) = (-17, -4.5, 3.1) for < 84 GeV and (-19, 0.68, -0.057) for > 84 GeV. 
This limit may be improved by an order of magnitude by the recently proposed GENIUS 
project ||122|| , or even by CRESST itself, assuming three years of operation with improved 
background rejection ||121| . 

Among high energy colliders, the Large Hadron Collider (LHC) at CERN is the ultimate 
supersymmetry discovery machine and will discover at least some superpartners in all of 
the regions of parameter space we have plotted. The LHC is scheduled to begin operation 
in its low luminosity mode in 2006. Before that, however, both the LEP II and Tevatron 
colliders have a chance to discover superpartners. The most stringent constraint from LEP 
II on minimal supergravity comes from chargino searches. LEP II is now concluding its run, 
and by the end of 2000 will improve the current chargino mass limits by about 5 GeV. If 
no signal is seen, this will marginally extend the bottom and right excluded regions of our 
figures. 

The Tevatron will begin operation early in 2001. In the first two years. Run Ha will 
provide an integrated luminosity of 2 fb~^ for each detector before a temporary shutdown 
for a year of detector maintenance and upgrades. In the subsequent Run lib, the data 
acquisition rate is expected to be about 5 fb~^ per year per detector. Hence, by 2006 we 
expect 10 — 12 fb~^ for each Tevatron collaboration. The Tevatron supersymmetry reach has 
been extensively studied recently [|123|| , with the conclusion that there is some sensitivity, but 
in rather limited regions of parameter space. The most effective signal is in the clean trilepton 
channel ||124| - |127[| resulting from chargino-neutralino pair production, followed closely by the 



jets plus channel ||128|| and the dileptons plus r jet channel ||129|| . The maximal reach 
in chargino mass is 170 GeV in the Bino LSP region at very low mo, where the leptonic 
branching ratios of the electroweak gauginos are enhanced by light sleptons. This degrades 
rapidly at higher mo, where hadronic decays are prominent. It also requires moderate tan/3. 
At large values of tan /?, decays to r leptons dominate the small mo region and signatures 
with r jets must be used 
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At present there are no dedicated Tevatron studies in the focus point region. (For an 
LHC study, see Ref. 



130 



There are several important modifications to coUider signals 
for mo > 1 TeV. For example, the lighter chargino and neutralinos are more degenerate, 
leading to softer decay products, and their branching ratios to b quarks are enhanced by 
their Higgsino component. Such issues may have a large impact on chargino and neutralino 
searches at the Tevatron. This is an important question, but currently the Tevatron reach 
in the focus point region is unknown. 

While the Higgs boson is not a supersymmetric particle, supersymmetry (in its economi- 
cal implementations) restricts its mass, and so Higgs boson searches also have an important 
impact on supersymmetric models. For LEP II, the ultimate exclusion limit, barring a dis- 
covery, is expected to be rrih > 115 GeV. At the Tevatron, the 3cr (5cr) Higgs boson discovery 
reach for 10 fb~^ is ruh ^ 100 (120) GeV [|131| - |133| . The Higgs boson mass, unlike all other 
quantities investigated here, is sensitive to the Aq parameter. As the Higgs boson mass limit 
rises, models with non-zero Aq, large tan/3, and mo ^ 1 TeV are increasingly favored ||22|| . 
However, for natural values of Aq p2|, 100 GeV < ruh < 120 GeV and so the Higgs boson 



will be discovered at the Tevatron at 3a, but never at 5a. 

Finally, there are many opportunities for discovering supersymmetry in low energy ex- 
periments. These include effects in hadronic and leptonic fiavor violation, CP violation, 
proton decay, and electric and magnetic dipole moments. These are discussed more com- 
pletely in Ref. ||2^. Here we will focus on two particularly robust probes: B Xgj and the 
anomalous magnetic moment of the muon. 

The best current measurements of B ^ Xs'~f from CLEO and ALEPH can be combined 
in a weighted average of B{B Xs'~f)exp = (3.14 ±0.48) x 10~^ ||134|| . These measurements 
will be improved at the B factories, where large samples of B mesons will greatly reduce 
statistical errors. However, the uncertainty in the theoretical prediction of the standard 
model, B{B X,7)sm = (3.29 ± 0.30) x 10-^ is likely to remain unchanged. By 2006, a 
conservative estimate is that both theoretical and experimental uncertainties will be ~ 0.3 x 
10"^. Combining them linearly, the 2a limit will be 2.1 x 10"^ < B{B Xsj) < 4.5 x lO"''. 

The supersymmetric contribution to the muon magnetic dipole moment (MDM) = 
^{g — 2)^ is also a robust probe, since it involves only a few (fiavor- and CP-conserving) 



|136|1 and is 



parameters [|135|] . The world average is a^^^P = (116 592 05 ± 45) x 10"^° 
consistent with the standard model. However, once data currently being taken is analyzed, 
the Brookhaven experiment E821 is expected to reduce the uncertainty to Aa^ ~ 4 x 
10~^° ||137|| . At present, uncertainties in the standard model prediction are substantial. 
Assuming these can be reduced, however, a reasonable estimate for future 2a sensitivity is 



In Table |^ we present our estimates for sensitivities that will be achieved before the 
LHC begins operation. The experiments likely to achieve these projections are also listed. 
Using these estimates, the reach in minimal supergravity parameter space for each mode 
is given in Figs. |l^ and In reading these figures, recall that we have assumed constant 
local densities in our assessment of dark matter search reaches. If one assumes that the 
local density is modulated by the thermal relic density, the dark matter reaches outside 
the shaded regions should be suitably diminished. Within the shaded regions, however, our 
analysis applies without modification. 

Several striking features emerge from Figs. IT8| and First, we see that, within the 
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TABLE IV. Constraints on super symmetric models used in Figs. 18 and 19 
experiments likely to reach these sensitivities before 2006. 



We also list 



Observable 



Type 



Bound 

> 100 GeV 

X, 



Experiment (s) 



X^X° 

Muon MDM 

''^proton 

V from Earth 
u from Sun 
7 (gal. center) 
7 (gal. center) 
cosmic rays 



Collider 
Collider 
Low energy 
Low energy 
Direct DM 
Indirect DM 
Indirect DM 
Indirect DM 
Indirect DM 
Indirect DM 



See Refs. il24|J127i^ 

\/\B{B X,7)| < 1.2 X lO-'' 



^ SUSY I 



< 8 X 10 



-10 



Equation ( pj|) 

< 100 km-2 yr-i 

< 100 km-2 yr-i 



yr 

(1) < 1.5 X 10-10 



cm 



^>-y(50) < 7 X 10-12 cm-2 s- 

(5/5)max < 0.01 



LEP: ALEPH, DELPHI, L3, OPAL 

Tevatron: CDF, DO 

BaBar, BELLE 

Brookhaven E821 

CDMS, CRESST, GENIUS 

AMANDA, NESTOR, ANTARES 

AMANDA, NESTOR, ANTARES 

GLAST 

MAGIC 

AMS-02 



minimal supergravity framework, nearly all of the cosmologically preferred models will be 
probed by at least one experiment. This is strictly true for tan/3 = 10. For tan/5 = 50, 
some of the preferred region escapes all probes, but this requires M1/2 ^ 450 GeV and 
mo ^ 1.5 TeV, and requires significant fine-tuning of the electroweak scale. In the most 
natural regions, all models in which neutralinos form a significant fraction of dark matter 
will yield some signal before the LHC begins operation. 

Also noteworthy is the complementarity of traditional particle physics searches and in- 
direct dark matter searches. Collider searches require, of course, light superpartners. High 
precision probes at low energy also require light superpartners, as the virtual effects of su- 
perpartners quickly decouple as they become heavy. Thus, the LEP and Tevatron reaches 
are confined to the lower left-hand corner, as are, to a lesser extent, the searches for devia- 
tions in i? — > Xs7 and a^. These bounds, and all others of this type, are easily satisfied in 
the focus point models with large mo, and indeed this is one of the virtues of these models. 
However, in the focus point models, all of the indirect searches are maximally sensitive, as 
the dark matter contains a significant Higgsino component. Direct dark matter probes share 
features with both traditional and indirect searches, and have sensitivity in both regions. It 
is only by combining all of these experiments, that the preferred region may be completely 
explored. 3 

Finally, these results have implications for future colliders. In the cosmologically pre- 
ferred regions of parameter space with 0.1 < VL^h"^ < 0.3, all models with charginos or 
sleptons lighter than 300 GeV will produce observable signals in at least one experiment. 



This is evident for tan/3 = 10 and 50 in Figs. O and 19. In Fig. 20, we vary tan/3, fix- 



ing Ml/2 to 400 GeV, which roughly corresponds to 300 GeV charginos. We see that the 



^Note that the complementarity referred to here is not the commonly recognized one, which 
concerns the mass of the neutralino. It is well-known that some indirect searches are effective even 
for LSP masses in the TeV range, well beyond the range of colliders. However, such models are 
highly unnatural, and they have not been considered here. 
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FIG. 18. Estimated reaches of various high-energy colhder and low-energy precision searches 
(black), direct dark matter searches (red), and indirect dark matter searches (blue) before the LHC 
begins operation, for tan /? = 10. The projected sensitivities used are given in Table 0. (The LEP 
chargino mass bound will marginally extend the bottom and right excluded regions and is omitted.) 
The shaded regions are as in Fig. |^. The regions probed extend the curves toward the forbidden, 
green regions. The dark matter reaches are not modulated by the thermal relic density. Bounds 
from photons from the galactic center are highly halo model-dependent; we assume a moderate 
halo profile parameter J = 500. (See text.) 



600 



500 



> 400 

cu 
O 



^ 300 



200 



100 



— 1 1 r- 



-I 1 r 








500 1000 1500 

mo (GeV) 



tan|8=50 



2000 



2500 



FIG. 19. As in Fig. 18, but for tan /? = 50. Here the <I>^ probe is sensitive to all of the parameter 
space shown and so its limit contour does not appear. 
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500 1000 1500 2000 2500 

mo (GeV) 

FIG. 20. As in Fig. but in the (mo, tan /3) plane for fixed M^/a = 400 GeV, = 0, and 
jj, > 0. The regions probed are toward the green regions, except for where it is between the 
two contours. The top excluded region is forbidden by limits on the CP-odd Higgs scalar mass. 

preferred region is probed for any choice of tan (3. (For extremely low tan/3 and mo, there 
appears to be a region that is not probed. However, this is excluded by current Higgs mass 
limits for Aq = 0. These limits might be evaded if Aq is also tuned to some extreme value, 
but in this case, top squark searches in Run H of the Tevatron will provide an additional 
constraint.) 

These results imply that if any superpartners are to be within reach of a 500 GeV lepton 
collider, some hint of supersymmetry must be seen before the LHC begins collecting data. 
This conclusion is independent of naturalness considerations. While our quantitative analysis 
is confined to minimal supergravity, we expect this result to be valid more generally. For 
moderate values of tan (3, if the dark matter is made up of neutralinos, they must either be 
light, Bino-like, or a gaugino-Higgsino mixture. If they are light, charginos will be discovered. 
If they are Bino-like, light sfermions are required to mediate their annihilation, and there 
will be anomalies in low energy precision measurements. And if they are a gaugino-Higgsino 
mixture, at least one indirect dark matter search will see a signal. For large tan (3, low energy 
probes become much more effective and again there is sensitivity to all probe superpartner 
spectra with light superpartners. Thus it appears, on qualitative grounds, that all models in 
which the scalar masses are not widely separated, and the charginos are not extravagantly 
heavy, will be accessible prior to LHC operation. 

VII. CONCLUSIONS 

In this paper, we have examined a wide variety of indirect dark matter detection signals. 
These include neutrinos from annihilation of dark matter in the cores of the Earth and Sun, 
continuum gamma rays from annihilation in the galactic center, and positron excesses in 
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cosmic rays from annihilation in the local solar neighborhood. In each case, the experimental 
landscape will be transformed in the next few years by experiments that are running or being 
mounted. We have tabulated the salient features and reaches of some of the most promising 
experiments in the previous sections. 

We have evaluated the prospects for dark matter detection in the framework of minimal 
supergravity. This framework incorporates many of the most compelling features of super- 
symmetry. Previously, this framework has been thought to predict a Bino-like LSP. That 
severely limited its utility for dark matter studies. However, recent work has made it clear 
that gaugino-Higgsino mixtures and even Higgsino-like LSPs are also quite naturally realized 
in minimal supergravity. We have been careful to include the full range of possibilities, with 
important (and positive) implications for future dark matter searches. 

Let us note in passing that in our parametrization of experimental probes, the case of no- 
scale supergravity ||138|| , recently revived in the context of gaugino-mediated supersymmetry 



breaking ||139| - |142|| , can be regarded as the special case tuq = 0. Experimental probes of these 
models are simply evaluated by restricting to the mg = axis. Several experiments, notably 
the trilepton Tevatron search, direct dark matter searches, and Brookhaven experiment E821 
will have the power to confirm or exclude this possibility in the near future. 

We have concentrated here on discovery signals. If a signal is confirmed, precision mea- 
surements may allow experiments to determine dark matter properties. For example, as 
has been noted in the literature, the neutralino's mass may be determined by the angular 
spread of the signal in neutrino telescopes. The energy spectrum of gamma rays or positrons 
signals may provide similar information. The gaugino-ness may also be constrained; indeed, 
the existence of a significant signal in itself would constitute evidence in favor of mixed 
gaugino- Higgsinos . 

The simplicity of minimal supergravity allows us to compare the reaches of a great variety 
of probes. We summarize by collecting several of our main conclusions: 

• Bino-like dark matter leads to suppressed rates for all indirect dark matter signals. In 
this case, unless the neutralino is extremely light (near current bounds), all indirect 
signals are beyond detection for the foreseeable future. 

• Higgsino-like dark matter cannot yield cosmologically interesting relic densities in a 
straightforward way. Studies that assume Higgsino-like dark matter exaggerate the 
power of indirect searches. 

• Mixed gaugino-Higgsino dark matter gives both relic densities in the preferred range 
0.1 ^ ^x^^ ~ 0-3 c'-iid detectable signals. Such dark matter is naturally present in 
focus point models, which are favored by low-energy constraints. 

• Experiments that are running or underway will transform the prospects for indirect 
dark matter detection. Among the most promising experiments are the neutrino 
telescopes AMANDA, NESTOR, and ANTARES; the gamma ray telescope MAGIC, 
and the satellite detector GLAST; and AMS-02, the anti-particle/anti- matter search 
aboard the International Space Station. For mixed gaugino-Higgsino dark matter, 
these experiments will be sensitive to nearly all models with cosmologically interest- 
ing neutralino relic densities, and are competitive with next-generation direct search 
experiments, such as CDMS (Soudan) and CRESST. 
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• The various indirect searches rely on different sources of neutralino annihilation (cores 
of the Earth or Sun, galactic center, local solar neighborhood) and so are sensitive to 
different assumptions. In addition, some signals, particularly the continuum photons, 
will be difficult to identify unambiguously as a dark matter signal. Without actual 
data and detailed analyses, it is difficult to make a more precise statement. However, 
we have seen that many experiments arc sensitive to the same supersymmetric models, 
and given the underlying uncertainties, redundancy is clearly a virtue. 

• Indirect searches are complementary to traditional particle searches. Separately they 
probe only portions of the cosmologically preferred model space. Combined, essentially 
all cosmologically preferred models will produce at least a hint of a signal in one of 
these experiments before the LHC begins operation. 

• In minimal supergravity models with 0.1 < fl^h'^ < 0.3, if there is no hint of supersym- 
metry before the LHC begins operation, no superpartners will be within reach of a 500 
GeV lepton collider. Our arguments are independent of naturalness considerations, 
and their qualitative structure suggests that similar conclusions will remain valid in 
alternative frameworks. 
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